Quantum computation with Josephson-qubits by using a current-biased information bus 
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We propose an effective scheme for manipulating quantum information stored in a superconducting nanocir- 
cuit. The Josephson qubits are coupled via their separate interactions with an information bus, a large current- 
biased Josephson junction treated as an oscillator with adjustable frequency. The bus is sequentially coupled to 
only one qubit at a time. Distant Josephson qubits without any direct interaction can be indirectly coupled with 
each other by independently interacting with the bus sequentially, via exciting/de-exciting vibrational quanta in 
the bus. This is a superconducting analog of the successful ion trap experiments on quantum computing. Our 
approach differs from previous schemes that simultaneously coupled two qubits to the bus, as opposed to their 
sequential coupling considered here. The significant quantum logic gates can be realized by using these tunable 
and selective couplings. The decoherence properties of the proposed quantum system are analyzed within the 
Bloch-Redfield formalism. Numerical estimations of certain important experimental parameters are provided. 

PACS numbers: 03.67.Lx, 74.50.+r, 85.25.Cp 



I. INTRODUCTION. 

The coherent manipulation of quantum states for realizing 
certain potential applications, e.g, quantum computation and 
quantum communication, is attracting considerable interest 
In principle, any two-state quantum system works as a 
qubit, the fundamental unit of quantum information. How- 
ever, only a few real physical systems have worked as qubits, 
because of requirements of a long coherent time and operabil- 
ity. Among various physical realizations, such as ions traps 
(see, e.g., OlHH]), QED cavities (see, e.g., S&l), quantum 
dots (see, e.g., |3 and NMR (see, e.g., j^, flOil ^ etc., su- 
perconductors with Josephson junctions offer one of the most 
promisin g pl atforms for reali zing qua ntum computation (see, 

l2al27ll2Ml29ll30ll3lll '). The nonlinearity of Josephson junc- 
tions can be used to produce controllable qubits. Also, cir- 
cuits with Josephson junctions combine the intrinsic coher- 
ence of the macroscopic quantum state and the possibility to 
control its quantum dynamics by using voltage and magnetic 
flux pulses. In addition, present-day technologies of integra- 
tion allow scaling to large and complex circuits. Recent ex- 
periments have demonstrated quantum coherent dynamics in 
the time domain in both single-qubit (see, e.g., jl2.fl3Lfl4ll ') 
and two-qubit Josephson systems L15J . 

There are two basic types of Josephson systems used to im- 
plement qubits: charge qubits |12] and flux qubits fTs"), de- 
pending on the ratio of two characteristic energies: the charg- 
ing energy Ec and the Josephson energy Ej. The charge 
qubit is a Cooper-pair box with a small Josephson coupling 
energy, Ej ^ Ec, and a well defined number of Cooper 
pairs is well defined. The flux qubit operates in another ex- 
treme limit, where Ej ^ Ec and the phase is well de- 
fined. A "quantronium" circuit operating in the intermedi- 
ate regime of the former two has also been proposed |14]. 
Voltage-biased superconducting quantum interference devices 



(SQUIDs), which work in the charge regime and with con- 
trollable Josephson energies, form the SQUID-based charge 
qubits that we will consider in this work. Our results can be 
extended to flux and flux-charge qubits. 

The key ingredient for computational speedup in quantum 
computation is entanglement, a property that does not ex- 
ist in classical physics. Thus, manipulating coupled qubits 
plays a central role in quantum information processing (QIP). 
Heisenberg-type qubit-couplings are common for the usual 
solid state QIP systems, e.g., the real spin states of the elec- 
trons in quantum dots |7, 8]. However, the interbit couplings 
for Josephson junctions involve Ising-type interactions, as su- 
perconducting qubits with two macroscopic quantum states 
provide pseudo-spin- 1/2 states. Recently, either the current- 
cuiTent interaction, by connecting to a common inductor, or 
the charge-charge coupling, via sharing a common capacitor, 
have been proposed to directly couple two Josephson charge 
qubits: the ith and jth ones. These interactions implement 
cri'^ crP^-type fisl fl^, ai"^ ® |17|, and the 

Ox ® cri -type 1 18] Ising couplings, respectively. Compared 
to the single-qubit operations, the two-qubit operations based 
on these second-order interactions are more sensitive to the 
environment. Thus, quantum decoherence can be more prob- 
lematic. In addition, capacitive coupling between qubits is 
not easily tunable 1 15|. Thus adjusting the physical parame- 
ters for realizing two-qubit operation is not easy. In order to 
ensure that the quanta of the relevant LC oscillator is not ex- 
cited during the desired quantum operations, the time scales 
of manipulation in the inductively coupled circuit should be 
much slower than the eigenfrequency of the iC-circuit 1 17]. 

Alternatively, the Josephson qubits may also be coupled 
together by sequentially interacting with a data bus, instead 
of simultaneously. This is similar to the techniques used for 
trapped ions ]2, 3], wherein the trapped ions are entangled by 
exciting and de-exciting quanta of their shared center-of-mass 
vibrational mode (i.e., the data bus). This scheme allows for 



faster two-qubit operations and possesses longer decoherent 
times. In fact, an externally connected LC-resonator (l^ and 
a cavity QED mode |201 were chosen as alternative data buses. 
However, it is not always easy to control all the physical prop- 
erties, such as the eigenfrequencies and decoherence, of these 
data buses. 

A large (e.g., up to lO/im) current-biased Josephson junc- 
tion (CBJJ) 1 2 1 1 is very suitable to act as information bus for 
coupling Josephson qubits. This because: i) the CBJJ is an 
easily fabricated device |22] and may provide more effective 
immunities to both charge and flux noise; ii) due to its large 
junction capacitance, the CBJJ can enable to be capacitively 
coupled over relatively long distances; iii) the quantum prop- 
erties, e.g., quantum transitions between the junction energy 
levels, of the current-biased Josephson junction are well estab- 
lished I23II24I1 ; and iv) its eigenfrequency can be controlled by 
adjusting the applied bias-current. In fact, a CBJJ itself can be 
an experimentally realizable qubit, as demonstrated by the re- 
cent observations of Rabi oscillations in them |25, 26J. Two 
logic states of such a qubit are encoded by the two lowest 
zero-voltage metastable quantum energy levels of the CBJJ. 
The decoherent properties of this CBJJ-qubit were discussed 
in detail in |27|. Experimentally, the entangled macroscopic 
quantum states in two CBJJ-qubits coupled by a capacitor 
were created |28|. Also, by numerical integration of the time- 
dependent Schrodinger equation, a full dynamical simulation 
of two-qubit quantum logic gates between two capacitively 
coupled CBJJ-qubits was given in |29]. 

In this paper, we propose a convenient scheme to selectively 
couple two Josephson charge-qubits. Here, a large CBJJ acts 
only as the information bus for transferring the quantum in- 
formation between the qubits. Thus, hereafter the CBJJ will 
not be a qubit, as in U^^^ ZL 2^ 2^- Two cho- 
sen distant SQUID-based charge qubits can be indirectly cou- 
pled by sequentially interacting these with the bus. This cou- 
pling method provides a repeatable way to generate entan- 
gled states, and thus can implement elementary quantum logic 
gates between arbitrarily selected qubits. Our proposal shares 
some features with the circuits proposed in 1 17, 18, 19, 21], 
but also has significant differences. Our proposal might be 
more amenable to experimental verification. 

The outline of the paper is as follows. In Sec. II we propose 
a superconducting nanocircuit with a CBJJ acting as the data 
bus, and investigate its elemental quantum dynamics. The bus 
is biased by a dc current and is assumed to interact with only 
one qubit at a time. There is no direct interaction between 
qubits. Therefore, the elemental operations in this circuit con- 
sist of: i) the free evolution of the single qubit, ii) the free evo- 
lution of the bus, and iii) the coherent dynamics for a single 
qubit coupled to the bus. In Sec. Ill we show how to real- 
ize the elemental logic gates in the proposed nanocircuit: the 
single-qubit rotations by properly switching on/off the applied 
gate voltage and external flux, and the two-qubit operations by 
letting them couple sequentially to the bus. The vibrational 
quanta of the bus is excited/absorded during the qubit-bus in- 
teractions. In Sec. IV we analyze the decoherent properties 
of the present qubit-bus interaction within the Bloch-Redfield 
formalism |32|, and give some numerical estimates for exper- 




Data Bus 



SQUID-Qubits 



FIG. 1 : SQUID-based charge qubits coupled via a large CBJJ. 



imental implementations. Conclusions and some discussions 
are given in Sec. V. 



II. A SUPERCONDUCTING NANOCIRCUIT AND ITS 
ELEMENTARY QUANTUM EVOLUTIONS. 

The circuit considered here is sketched in Fig. 1. It con- 
sists of N voltage-biased SQUIDs connected to a large CBJJ. 
The fcth (fc = 1, 2, N) qubit consists of a gate electrode 
of capacitance Cg^ and a single-Cooper-pair box with two ul- 
trasmall Josephson junctions of capacitance Cj^ and Joseph- 
son energy E'j^ , forming a DC-SQUID ring. The inductances 
of these DC-SQUID rings are assumed to be very small and 
can be neglected. The SQUIDs work in the charge regime 
with fcfiT <C Ej ^ Ec ^ A, in order to suppress quasi- 
particle tunneling or excitation. Here, fc^. A, Ec, T, and Ej 
are the Boltzmann constant, the superconducting gap, charg- 
ing energy, temperature, and the Josephson coupling energy, 
respectively. 

The connected large CBJJ biased by a dc current works in 
the phase regime with Ej ^ Ec- It acts as a tunable an- 
harmonic LC-resonator with a nonuniform level spacing and 
works as a data bus for transferring quantum information be- 
tween the chosen qubits. The mechanism for manipulating 
quantum information in the present approach is different from 
that in I2TI1 . although the circuit proposed here 

might seem similar to those there. The differences are: 

(1) a large CBJJ, instead of LC-osciflator O El El 
formed by the externally connected inductance L and the ca- 
pacitances in circuit, works as the data bus; 

(2) we modulate the applied external flux, instead of the 
bias-current |21|, to realize the perfect coupling/decoupling 
between the chosen qubit and the bus; and especially 

(3) the free evolution of the bus during the operational de- 
lays will be utilized for the first time to control the dynamical 
phases for implementing the expected quantum gates. 

The Hamiltonian for the present circuit can be written as 
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Here, n^, = C,,Vk/i2e), Ck = C,, + Cj,, Cj, = 20%, 

Cb = Cb + Eti Cj,C,JCk,Ej, = 2El cos(7r<i>,./<i>o), 
and 9k = {Ok2 + ^fci)/2 with Ok^ and being the phase 
drops across two small Josephson junctions in the fcth qubit, 
respectively. Also, Cg^, $o> ^k, and Vk are the gate capac- 
itance, flux quantum, external flux, and gate voltage applied 
to the fcth qubit, respectively. Correspondingly, Cf,, 0b, Eb, 
and lb are the capacitance, phase drops, Josephson energy, 
and the bias-current of the large CBJJ, respectively. Above, 
the number operator nj. of excess Cooper-pair charges in the 
superconducting island and the phase operator 9k of the or- 
der parameter of the fcth charge qubit are a pair of canonical 
variables and satisfy the commutation relation: 

[9k,nk] = i. 

The operators db and p f, are another pair of canonical variables 
and satisfy the commutation relation: 
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with 2TTpb/^o = 2nbe representing the charge difference 
across the CBJJ. 

The CBJJ works in the phase regime. Thus, Ec,, = 
e'^/{2Cb) <C Eb and the quantum motion ruled by the 
Hamiltonian Hr equals to that of a particle with mass m = 
C'b($o/27r)2 in a potential U{9b) = -Eb{cos9b + hOb/Ir), 
Ir = 2TrEb/^o. For the biased case h < Ir, there ex- 
ists a series of minima of U{9b), where dU{9b) / d9b = 
0, d'^U{9b)/d9l > 0. Near these points, U{9b) approximates 
to a harmonic oscillator potential with a characteristic fre- 
quency 
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depending on the applied bias-current Correspondingly, 
the Hamiltonian Hr reduces to 
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number of quantum metastable 
of the quantum oscillator is 

1 - Ib/Irf^^ . 



The energy scale of the quantum oscillator (3) is LUb / (27r) ^ 
10 GHz [25], which is of the same order of the Josephson 
energy in the SQUID. Therefore, the oscillating quantum of 
the information bus will be really excited, even if only one of 
the qubits is operated quantum mechanically. This is different 
from the case considered in wherein the LC-oscillator 
shared by all charge qubits are not really excited, as the eigen- 
frequency of the LC-circuit is much higher than the typical 
frequencies of the qubits dynamics. For operational conve- 
nience, we assume that the bus is coupled to only one qubit 
at a time. The coupling between any one of the qubits (e.g., 
the fcth one) and the bus can, in principle, be controlled by 
adjusting the applied external flux (e.g., In this case, any 
direct interaction does not exist between the qubits, and the 
dynamics of the CBJJ can be safely restricted to the Hilbert 
space spanned by the two Fock states: |0b) and \lb), which 
are the lowest two energy eigenstates of the harmonic oscil- 
lator of Eq. (3). Furthermore, we assume that the applied 
gate voltage of any chosen (fcth) qubit works near its degen- 
eracy point with Ug^ = 1/2, and thus only two charge states: 
\nk = 0) = I tfc) and \nk = 1) = | ik), play a role during 
the quantum operation. All other charge states with a higher 
energies can be safely ignored. Therefore, the Hamiltonian 
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describes the interaction between any one of the qubits (e.g., 
the fcth one) and the bus, and provides the basic dynamics for 
the present network. Here, SEc^ — 2e^(l — 27ig^)/Cfc, Afc — 

Cg^ (2n/^o)^h/{2CbLOb)/{2Ck), and the pseudospin opera- 
tors are defined by: 



IU)(UI + Ufe)(Tfe|, 
-*ITfe)(UI + *Ufe)(Tfc| 
\U){U\-\ik){lk\. 



Above, when the first cosine-term in Hamiltonian (2) was ex- 
panded, only the single-quantum transition process approxi- 
mated to the first-order of 9b was considered. The higher or- 
der nonlinearities have been neglected as their effects are very 
weak. In fac t, for the lower number states of the bus, we have 

eters jllf ll^ 

Notice that the coupling strength Afc between the qubit and 
the bus is tunable by controlling the flux $fc, applied to the 
selected qubit, and the bias-current /;,, applied to the informa- 
tion bus. For example, such a coupling can be simply turn off 
by setting the flux "I>fc as $o/2. This allows various elemen- 
tal operations for quantum manipulations to be realizable in a 



2, for the typical experimental param- 

Cb ^ IpF, Wb/27r - 10 GHz, and 



4 



controllable way. In the logic basis {|Ofe), | Ife)}, defined by 

^ I ifc) + I Tfc) ^ I ife) - I Tfc) 



v/2 



and under the usual rotating-wave approximation, the above 
Hamiltonian (4) can be rewritten as 
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with 



|U)(0fe| + |0fe)(U|, 
-i|lfc)(Ofc|+z|Ofe)(lfc| 



and (j'^^ = {ai''^ ± (Ty'^-')/2. Here, the logic states \0k) and 
1 1/c) correspond to the clockwise and anticlockwise persistent 
circulating currents in the fcth SQUID-loop, respectively. 

We now discuss the quantum dynamics of the above 
Josephson network. Without loss of generality, we assume 
in what follows that the bias-current lb applied to the CBJJ 
doesn't change, once it is set up properly beforehand. The 
quantum evolutions of the system are then controlled by other 
external parameters: the fluxes applied to the qubits and the 
voltages across the gate capacitances of the qubits. Depending 
on the different settings of the controllable external parame- 
ters, different Hamiltonians can be induced from Eq. (6) and 
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thus different time-evolutions are obtained. Obviously, during 
any operational delay r with c&jf. = $o/2 and Vk = e/Cg., 
the ith qubit remains in its idle state because the Hamiltonian 
vanishes (i.e., H^^ = 0) as Ej^ = 0, Ug. = 0. However, the 
data bus still undergoes a free time-evolution: 



Uo{t) = exp ( -^Hb 



(7) 



This evolution is useful for controlhng the dynamical phase of 
the qubits to exactly realize certain quantum operations. For 
the other cases, the dynamical evolutions of the chosen qubit 
depend on the different settings of the experimental parame- 
ters. 

1) For the case where = <l^(}/2 and Vk ^ e/Cg,^, the 
«th qubit and the bus separately evolve with the Hamiltonians 
h['''^ = -5Ec^df^ /2 and determined by Eq. (3), re- 
spectively. The relevant time-evolution operator of the whole 
system reads 



C/f^t) = exp 



) exp 



(8) 



2) If the /cth qubit works at its degenerate point and couples 
to the bus, i.e., Vk = e/Cg^ and $fc ^ $o/2, then we have 
the Hamiltonian 



Hkb = Ej,af^/2 + Hb + iXk aa^+^ 



a^ai'^ (9) 



from (6). The corresponding dynamical evolutions are 



|0fa)|0fe) ^ e*^'=*/2|0fa)|0,), Ukb = eM-iHkbt), = Ejjh-cob, 
\Ob)\lk) ^ e--"*! [cos(%t) -i^sin(%t)] |0,)|lfe) -^sm{^t) \h)\Ok)} 
^ \hm) ^ e--"* { [cos {^t)+i^sm{^t)] |l,)|Ofe) + |^sin(^^i) |06)|1,)} 
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wifli Qk = V^l + (2Afe/fi)2. 



Specifically, we have the time-evolution operator 
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with 



A{t) = exp 




for the resonant case: = 0. This reduces Eq. (10) to the 
time evolutions: 
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\0b)\0k) |06)|0fe), 

\Ob)\W ^Ai.e--'.* [cos(Y) |0,)|U)-sin(Y)|l6>|0.)], 
[ IWm^^e-^"' [cos(Y)|l6)|0fe)+sin(^)|0fe)|lfe)]. 



For another extreme case, i.e., the system works in the dis- with 
persive regime (far from the resonant point): 2Afc /{h\Ak\) <C 

1, we have the time evolution operator £r/ ,2/11 \/i i--t in \/n l-t~^//fcA \ 

U^''^ {t) = A{t) exp -^1^^ . (12) It reduces to the following time evolutions: 
I 



\0b)\0k) exp(zf^)|0(,)|0,), 



\Wk} 
\h)\Ok) 
\h)\h) 



exp 



exp 



exp 



3) Generally, if $fe ^ $o/2 and Fg^, ^ e/Cg^, then the 
Hamiltonian (6) can be rewritten as 



4k) ^ 
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and (T^'^-' = (o-i*^^ ± (j''y^')/2. Here, 003% = Ej^/E]^, and 
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$fe are set properly beforehand such that Ej,^ ~ HuJt, <^ 
SEcf, , then the detuning HAk = Ek~ hujb is very large (com- 
pared to the coupling strength Afe < 10~^Eji^). Therefore, the 
time-evolution operator of the system can be approximated as 



^i*=Ht)=S(Oexp|-i 



with 
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^{SEci^y + Ej^. If the bias-current lb and the flux This implies the foUowing evolutions 
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\0b)\0k) e-<'=*{[cos(aO + »cos77fcSin(^fci)] |Ofa) |Ofc) + i sinryfe sin(^feO |Ob)|lfc)}, 

|06)|lfc) ^^^^ {[cos(efei) -z cos 77feSin(Cfei)] |Ofc) |lfc) + i sin?7fe sin(Cfci) |06)|0fc)} , 

|l6)|0fc) e-'(«'= {[cos(^;,<) +icosr,k sin(e^t)] |l6)|0fc) + zsinr;^ sm(e^t) |lfc)|lfe)} , 

. e-'(f^+-'')*{*sin%sin(^^t) |lfc)|Ofc) + [cos(e^i) - z cost?, sin(^^t)] \h)\lk)} , 



with 
and 

ei. = a- + A^/(?i'A,). 

In what follows we shall show that any process for ma- 
nipulating the quantum information stored in the present cir- 
cuit can be effectively implemented by selectively using the 
above elementary time-evolutions: Uo{t), u[''\t), ij^\t), 

III. QUANTUM MANIPULATIONS OF THE 
SUPERCONDUCTING NANOCIRCUIT. 

It is well known that any valid quantum transformation can 
be decomposed into a sequence of elementary one- and two- 
qubit quantum gates. The set of these gates is universal, and 
any quantum computing circuit comprises only gates from this 
set. Several schemes 1171 [1^12^ have been proposed for im- 
plementing one of the universal two-qubit gates with Joseph- 
son qubits by using the direct interactions between them. By 
making use of the data bus interacting sequentially with the se- 
lective qubits, Blais et al. [21] showed that the two-qubit gate 
may be effectively realized. Two important problems will be 
solved in our indirect-coupling approach: 

i) when one of two qubits is selected to couple with the data 
bus, how we can let the remainder qubit decouple completely 
from the bus; and 

ii) the phase changes of the bus' and qubit's states during 
the operations are very complicated, how we can control these 
phase changes in order to precisely implement the desired 
quantum gate. 

The scheme in i2lll assumed that, when one of the two 
qubits is tuned to resonance with the bus, then the other qubit 
is hardly affected because of its different Rabi frequency. Ob- 
viously, this decoupling is not complete and thus it is not easy 
to assure that the bus couples only one qubit at a time. By 
controlling the external flux $fe applied to the qubits, the net- 
work proposed here provides an effective method for making 
the remainder qubit completely decouple from the bus. All the 
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desired elementary operations for quantum computing can be 
exactly implemented by properly setting the experimentally 
controllable parameters, e.g., the external the gate volt- 
age Vk, the bias-current /{,, and the duration t of each selected 
quantum evolution, etc. 

Hereafter, we assume that each of the selected time- 
evolutions can be switched on/off very quickly. 

A. single-qubit operations 

First, we show how to realize the single-qubit operations 
on each SQUID-qubit. This will be achieved by simply turn- 
ing on/off the relevant experimentally controllable parame- 
ters. For example, if Ug^ ^ 1/2 and Ej^ = for a time 

span t, then the time-evolution C/}*^' (t) in equation (8) is real- 
ized. This operation is the single-qubit rotation around the x 
axis: 

(cos ^ i sin ^ \ 
(15) 
i sin ^ cos ^ / 

with (pk = SEc^.t/h. Rotations by ifi = vr and ipk = 7''/2 
produce a spin flip (i.e., a NOT-gate operation) and an equal- 
weight superposition of logic states, respectively. 

The rotation around the z axis can be implemented by using 
the evolution (12). This operation is conditional and depen- 
dent on the state of the bus. If the bus is in the ground state 
|Ofc), the rotation reads 

i?W(0.) = e-^'^'*(^""^^' (16) 

with Qk = Wfc/2 + Xl/{2h'^Ak), 4>k = EjJ/{2h) + 
X\t / {21t? Ak) ■ With a sequence of x- and z-rotations, any 
rotation on the single-qubit can be performed. For example, 
the Hadamard gate applied to the fcth qubit: 

= 7i ( 1 -1 ) ' 

can be implemented by a three-step rotation: 

(17) 
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Here, the relevant durations ti, t2, and are set properly to 
satisfy the conditions 
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B. two-qubit operations 

Second, we show how to realize two-qubit gates by letting 
a pair of qubits (the fcth- and jth ones) interact separately with 
the bus. Before the quantum operation, the chosen qubits de- 
couple from the bus. At the end of the desired gate opera- 
tion the bus should be disentangled again from the qubits, and 
returned to its ground state. For operational simplicity, we 
assume that the bus resonates with the control qubit, the fcth 
one, i.e., = 0. We now consider the following three-step 
operational process: 

i) Couple the control qubit to the bus (i.e., the applied exter- 
nal flux is varied to $o) and realize the evolution u[''^ (ti) 
for the duration ti : 



sm 



-1. 



(18) 



Then, by returning the $fc to its initial value, i.e., = $o/2, 
the fcth qubit can be decoupled from the bus exactly. Before 
the next step operation, there is an operational delay ti . Dur- 
ing this delay the state of the qubits does not evolve, while the 
data bus still undergoes a time-evolution C/o(ti). 

ii) Couple the target qubit (the jth one) to the bus and real- 
ize the time-evolution J/g''^ {t2 ) . This is achieved by letting the 
chosen qubit work near its degenerate point (i.e., Ug. ^ 1/2) 
and switching on its Josephson energy (i.e., $j 4'o/2). Af- 
ter the time t2 determined by the condition 



C0s{^jt2) = - Sm{^jt2) = 1, 



(19) 



we decouple the jth qubit from the bus and let it be in the 
idle state by returning its gate-voltage Vj to the degenerate 
point (rig^. = 1/2), and simultaneously switching off the rel- 
evant Josephson energy. During another operational delay T2 
before the next step operation, the bus undergoes another free- 
evolution Uq{t2). 

iii) Repeat the first step and realize the evolution C/}'^' {ts) 
with 



sm 



1. 



(20) 



DiagranmiaticaUy, the above three-step operational process 
with two delays can be represented as foUows: 
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ie-'^-''''^{cos7]j\QblkOj) + smr]j\Oblklj)), 



le 

r'rW 



iX-ia;.(ti+t2+ri+r2)(gij^^^.|1^0feO,) - COS r/, |UOfclj)) 

ze-'>^-«^''^(sinr;, |06lfc0,) - cosr/. lOfclfclj}), 
I 



with T = ti + t2 + + Ti + T2 being the total duration operations. If the total duration T is satisfied as 
of the process, and x = Cjt^ + Wf,(Ti -|- r2)/2. Obviously, 

the information bus remains in its ground state |0(,) after the sin(a;6T) = 1, 



(21) 



the above three-step process with two delays yields a two- 
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qubit gate expressed by the following matrix form 



1 

1 

cosrij 

sin 




(22) 



which is a universal two-qubit Deutsch gate (3. 

Analogously, if the second step operation tj^^\t2) in the 

above three-step process is replaced by the operation 1)2^ (<2), 
then another two-qubit operation expressed by 



with Tj = exp{kjt2),Aj = exp(i^jt2), ^^j = Ej-/{2h) + 
X'j/{2h'^Aj), <rj = + \'jt2/{fi'^Aj), can be implemented. 
This three-step operational process can similarly be repre- 
sented diagrammatically as 



&^'''\t2) = 






V 





j 












(23) 



lOfcOfeOj 



\ObOklj) 



loafeO, 



lOhlfelj) 



C/o(ri)C/f' (ti) 



i>o(ri)C/f' (ti) 



T*e- 



lObOfcO,) 



C/f (t3) 



re--|0f,0fc0j), 



-»u.,(ti+3ri/2)|^^Q^Q^.\ 



C/o(r2)t/<^'(i2) 



C/f (t3) 



C/o(r2)t/«'(t2) 



il/-icJ6(tl+t2+Ti+T2) I IjjOfelj) 



with = LUbt2/2 + \jt2l{2h^Ai) + W6(ti + T2)/2. Above, 
the durations of the first- and third-step operations have been 
set the same as those for realizing the two-qubit operation 

The two-qubit gate uf''\T]j) (or U^^\t2)) performed 
above forms a universal set. Any quantum manipulation can 
be implemented by using one of them, accompanied by ar- 
bitrary rotations of single qubits. Obviously, if the system 
works in the strong charge regime: Ej./{SEci) ^ 1, and 

cos77j ~ 0, sin r/j ^ 1, then the two-qubit gate u'^''\ilj) in 
(22) approximates the well-known controUed-NOT (CNOT) 
gate 



known controUed-phase (CROT) gate 



CNOT ~ 



/ 1 
10 
1 

Vo 1 



Also, if the duration ^2 of the evolution J/2"'' (^2) and the de- 
lays Ti , T2 are further set properly such that 

cos(<;j t2) = sin(?jt2) = sin(tJbT) = 1, 
then the two-qubit operation U2 in (23) reduces to the well- 



U, 



CROT 



/ 1 
10 
10 

Vo -1 



IV. DECOHERENCE OF THE QUBIT-BUS SYSTEM DUE 
TO THE BIASED VOLTAGE- AND CURRENT-NOISES 

An ideal quantum system preserves quantum coherence, 
i.e., its time evolution is determined by deterministic re- 
versible unitary transformations. Quantum computation re- 
quires a long phase coherent time-evolution. In practice, any 
physical quantum system is subject to various disturbing fac- 
tors which destroy phase coherence. In fact, solid-state sys- 
tems are very sensitive to decoherence, as they contain a 
macroscopic number of degrees of freedom and interact with 
the environment. However, coherent quantum manipulations 
of the qubits are still possible if the decoherence time is fi- 
nite but not too short. Hence, it is important to investigate 
the effects of the environmental noise on the present quantum 
circuit. 
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II 





FIG. 2: Schematic diagram of a SQUID-based charge qubit with 
impedance Z{u)i) coupled to a CBJJ with admittance Y{uj2). 



The typical noise sources in Josephson circuits consist of 
the Hnear fluctuations of the electromagnetic environments 
(e.g., circuitry and radiation noises) and the low-frequency 
noise due to fluctuations in various charge/current channels 
(e.g., the "background charge" and "critical current" ). Usu- 
ally, the former one behaves as Ohmic dissipation [35] and the 
latter one produces a 1// spectrum js^. Within the present 
work, we will consider the case of Ohmic dissipation due to 
linear fluctuations of the external circuit parameters: the bias- 
cuiTent lb applied to the CBJJ and the gate voltages applied 
to the qubits. The effect of gate-voltage noise on a single 
charge qubit and that of bias-current noise on a single CBJJ 
has been discussed in 35] and in It27i] . respectively. We 
now study these noises together (see figure 2), since the in- 
teraction between a CBJJ, acting as a bus here, and a selected 
(e.g., the fcth) qubit takes a central role in the present scheme 
for quantum manipulations. Each electromagnetic environ- 
ment is treated as a quantum system with many degrees of 
freedom and modeled by a bath of harmonic oscillators. Fur- 
thermore, each of these oscillators is assumed to be weakly 
coupled to the chosen system. The Hamiltonian of a chosen 
(fcth) qubit coupling to the bus, containing the fluctuations of 
the applied gate voltage Vk and bias-current /;,, can be gener- 
ally written as 



H = Hkb+HB + V 



with 



- EE 



2m,. 



E E + J 



and 

V ^ 



sinafccr^'"'-' + cos akcr^^^ 



2 

huj-i 



(24) 



(25) 



being the Hamiltonians of the two baths and their interactions 
with the non-dissipative qubit-bus system H^b, respectively. 



Above, 



are the Boson operators of the jth bath, and 



Ri 



Ck 



i?2 



2CbL0b 



with g^. being the coupling strength between the oscillator 
of frequency u)j and the non-dissipative system. The effects 
of these noises can be characterized by their power spectra, 
which in turn depend on the corresponding "impedance" (or 
"inductance") and the temperature of the relevant circuits. For 
example, introducing the impedance Zt{uj) = l/[iujCt + 
Z^^{uj)] with Z{uj) — Rv being the Ohmic resistor, the 
corresponding voltage between the terminals of impedance 
Zt{ijj) can be expressed as 6V — A^^jX^j. Thus, the 
spectral density of this voltage source for Ohmic dissipation 
can be expressed as 



•E 



5{u) — uii) 



(26) 



Similarly, the spectral density for the bias-cuiTent source can 
be approximated as 



F{u) 



(27) 



with Yf being the dissipative part of the admittance of the cur- 
rent bias. 



The well-established Bolch-Redfield formalism I32ll37r of- 
fers a systematic way to obtain a generalized master equa- 
tion for the reduced density matrix of the system, weakly in- 
fluenced by dissipative environments. A subtle Markov ap- 
proximation is also made in this theory such that the resulting 
master equation is local in time. Of course, in the regime of 
weak bath coupling and low temperatures, this theory is nu- 
merically equivalent to a full non-Markovian path-integral ap- 
proach 1 38]. For the present qubit-bus system and in the basis 
spanned by the eigenstates {\g), |un), |w„),n — 1,2, ...} of 
the non-dissipative Hamiltonian Hkb, the Bloch-Redfield the- 
ory leads to the following master equations 



dt 



+ 2^ (^Q/3^iy + •^Q/J^iiy) cr^i/, (28) 



with 



10 



and 



R. 



1 f°° I 

^ _ V K 



(29) 



(30) 



with 



' g,{±T) = (^)'e.. l5..P[(n(a;i) + l)e^-- + (n(a;i))e±-^-], 
^ ^2(±r) = (^) E., l5u..P(n(a;2) + l)e^-^^ 



Above, each one of the states |a), ... can be equal to one 

of the eigenstates of Hkb- = l/[exp{HLUj/kBT) — 

1] is the average number of thermal photons in the mode of 
frequency Wj. The denotation = (a|x|/J) accounts for the 
matrix element of operator x, i.e., 

Ac,0 = {a\Ak\P), ife = ai'^hinak+ai^^ cosak = <jf\ 
and 

B^p = (a|a|/3), B^^ = {a\a^\p). 

Also, u>ai3 = [Ea — Eff) /H with Ea {Ejs) being one of eigen- 
values of the non-dissipative Hamiltonian Hkb, correspond- 
ing to the eigenstate \a) ( |/3) ). The spectrum of Hkb includes 



the ground state \g) = 0), corresponding to the energy 
Eg = —hAk/2, and a series of dressed doubled states 

\un) = cos6n\+k,n) -ism6n\-k,n+l), 

\vn) = -isin^„|-|-fc,n) -|-cos^„|-fe,n-|-l) 

corresponding to the eigenvalues 

Eun = t^bin + 1) - E^n = ^h{n + 1) + y , 
with 

cos6'„ = pn- hKkl^J {p„ - KKkY + 4Xlin + 1), 
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and 



^(?lA02+4A|(n + l). 



Here, |±fc) and |n) are the eigenstates of the operators ai'^^ 
and Hb with eigenvalues ±1 and hujb{n + 1/2), respectively. 

Under the secular approximation, the evolution of the non- 
diagonal element aaj3 of the reduced density matrix a is de- 
termined by 



di 



0-al3 + {i + Im(i?a/3a/3) + Im(S'a/3a/3)] 

+ [Re{Ra0a0) + ReCS'c/Ja/s)]} O-a/3 = 0. (31) 



Here, RaPfi,^ and Safifi,^ are calculated respectively from 
Rai^nv and Sajdixv by setting /x = a and v = (}. Re(a;) and 
Im(x) represent the real- and imaginary parts of the complex 
number x. The formal solution of the above differential equa- 
tion (31) reads 

o-a/3(i) = o-a/3(0) exp {-T'^t^ exp {-iQa^t) , (32) 



with Qap = ujap + lm{B.a(iai:i) + lm{Sapap) being the ef- 
fective oscillating frequency (the original Bohr frequency 
plus the Lamb shift Auaff = ImRa^a^ + ImSa^aff), and 



= -[R'e(-Ra/3a/3) + 'Re{Sa0a0)] 



(33) 



describing the rate of decoherence between the states \a) and 



In the present qubit-bus system operating near the resonant 
point: Ek ^ /iwh, the decoherences relating to the lowest three 
energy eigenstates, i.e., jg), jug) — |u),and j^o) = \v), are 
specially important for the desired quantum manipulations. 
The decoherences outside these three states are negligible. 
After a long but direct derivation, we obtain the decoherence 
rates of interest: 



av |4 (sinafc cos^6lo)^ ~~f~" ^ (cosafe cos6lo)^ coth ^ ^"^ ^ w„g 
-I- (cosafesin^o)^ coth ^ ^"^ ^ ~ 1 w^g + (sinafe sin2^o)^ coth ^ ^"y ^ ~ 1 



+ (cos a cos ^o) 



^ , 2 (cos a sin Oof coth ^ ^^^ ^ t^^,, 



+ a/ 



-|- (sin a sin 2^0 )^ 



coth 



.2A;bT, 



vgi 



IjJv 



(34) 



(35) 



and 



ay •{ 4 (sin a cos 26*0)^ — j h 2 (sin a sin 26*0)^ coth 



-|- (cos a cos ^0) 



coth 



h 

2kBT 



WbT) 



+ 1 



^^ug + (cos a sin ^o) 



coth 



+ aj < sin ^0 



coth ( + 1 



Ulug + COS 9o 



COthf ^'^^ 



2kBT I 







( f^vg \ 




bTJ 


+ 1 


Ulyg^ 



(36) 



Above, the various Bohr frequencies read 

: LOb/2 + Eu/{2h) - ^J{hiOb-Ek)^+4Xl/{2h), 



"vg 



tOb/2 + Ek/{2h) + ^{nwb-EkY+A\l/{2h), 



Jug 
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and 

w„„ = ^ {hu!b - EkY + 4:Xl/h. 

Two dimensionless parameters ay = 
ttRvC^J[RkCI], Rk = h/e^ w 25.8 kfl and 
a J = Yi/{CbiOb) characterize the coupling strengths 
between the environments and the system. 

Specially, if the system works far from the resonant point 
(with Afc ~ 0, achieved by switching off the Josephson en- 
ergy), the above results (shown in Eqs. (34-36)) reduce to 
those for the case when the qubit and the bus in- 

dependently decohere. Namely, reduces to the rate II 111 

which describes the decoherence between two charge states 
I i) and I t) of the superconducting box with zero Josephson 
energy. Also, T^^^ reduces to the decoherent rate 12711 

To,^ = ai[cothir2UJb/2kBT) + l]ujb, 

between the ground and first excited states of the data 
bus. However, for the strongest coupling case (i.e., when 
the system works at the resonant point), we have Ek = 
Ej^ = huJb,cosak — 1, cos 00 — sinOo = l/\/2, and 
coth[nuJug/l2kBT)] - 1 ~ coth[ftw„g/(2fcsT)] - 1 ~ 
(< 10~^, for the typical experimental parameters fT2il : 
Afe ~ O.lEj,, Ej, = hujb bOfieV > keT ~ S^eV). 
Thus, the minimum decoherent rates 

Tg'u = (av + ai)ujug, (37) 

Tgv^ = {av + ai)ujvg, (38) 

and 

uv gu ' gv ^ ^ / 

are obtained for the above three dressed states, respectively. 

It has been estimated in Ref. 1 11 ] that the dissipation for 
a single SQUlD-qubit is sufficiently weak: ay ^ 10^^ for 
Rv = 5Qfl, Cj^/Cgf. ^ 10~^, which allows, in principle, for 
10^ coherent single-qubit manipulations. For a single CBJJ 
the dimensionless parameter aj only reaches lO^'^ for typical 
experimental parameters f2^: l/Yj ^ 100 Q, Cb ^ 6 pF, 
LObl^TT ^ 10 GHz. This implies that the quantum coherence 
of the present qubit-bus system is mainly limited by the bias- 
current fluctuations. Fortunately, the impedance of the above 
CBJJ can be engineered |25| to be 1/Y/ ~ 560 kil. This lets 
aj reach up to 10^^ and allow about 10^ coherent manipula- 
tions of the qubit-bus system. 



V. CONCLUSIONS AND DISCUSSIONS 

In summary, we have proposed an effective scheme to cou- 
ple any pair of selective Josephson charge qubits by letting 



them sequentially couple to a common CBJJ, which can be 
treated as an oscillator with adjustable frequency. Two logic 
states of the present qubit are encoded by the clockwise and 
anti-clockwise persistent circuiting currents in the dc SQUID- 
loop. At most one qubit can be set to interact with the bus 
at any moment. The interaction between the selected qubit 
and the data bus is tunable by controlling the flux applied to 
the qubit and the bias-current applied to the data bus. This 
selective coupling provides a simple way to manipulate the 
quantum information stored in the connected SQUlD-qubits. 
Indeed, any pair of selective qubits without any direct interac- 
tion can be entangled by using a three-step coupling process. 
Furthermore, if the total duration is set up properly, the de- 
sired two-qubit universal gates, which are very similar to the 
CNOT- and CROT gates, can be implemented via such three- 
step operational processes. During this operation, the mode of 
the data bus is unchanged, although its vibrational quantum is 
really excited/absorbed. After the desired quantum operation 
is performed on the chosen qubits, the data bus disentangles 
from the qubits and returns to its ground state. 

In previous schemes, the distant Josephson qubits are cou- 
pled directly by either the charge-charge interaction, via con- 
necting to a common capacitor, or by a current-current in- 
teraction, via sharing a common inductor The present indi- 
rect coupling scheme offers some advantages: i) the coupling 
strength is tunable and thus easy to be controlled for realizing 
the desired quantum gate, ii) this first-order interaction is more 
insensitive to the environment, and thus possesses a longer 
decoherence time. Also, compared to pre vious data buses, 
the externally connected LC-resonator I19ll and cavity QED 
mode 1 20], the present CBJJ bus might be easier to control for 
coupling the chosen qubit. For example, its eigenfrequency 
can be controlled by adjusting the applied dc bias-current. In 
addition, the CBJJ is easy to fabricate using current technol- 
ogy i22il and may provide more effective immunities to both 
charge and flux noise. 

By considering the decoherence due to the linear fluctu- 
ations of the applied voltage T4 and current lb, we have 
analyzed the experimental possibility of the present scheme 
within the Bloch-Redfield formalism. A simple numeri- 
cal estimate showed that the quantum manipulations of the 
present qubit-bus system are experimentally possible, once 
the impedance Yi of the CBJJ can be engineered to have a 
sufficient low value, i.e., can be enlarged sufficiently 

(e.g., 1/Y/ - 560 K ^2 Q). Of course, this possibility, 
like those in previous schemes 113 El d US l21ll . is also 
limited by other technological difficulties, e.g., suppress the 
low-frequency 1// noise, and fast switch on/off the external 
flux to couple/decouple the chosen qubit, etc.. For example, 
a very high sweep rate of magnetic pulse (e.g., up to ^ 10^ 
Oe/s 1 39]), is required to change half of flux quantum through 
a SQUlD-loop (with the size e.g., 50/im) in a sufficiently short 
time (e.g., the desired ~ 40 ps). This and other obstacles pose 
a challenge that motivate the exploration of novel circuit de- 
signs that might minimize some of the problems that he ahead 
in the future. 
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